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Abstract 

For a disordered two-dimensional model of a topological insulator (such as a Kane-Mele 
model with disordered potential) with small coupling of spin invariance breaking term (such 
as the Rashba coupling), it is proved that the spin edge currents persist provided there is a 
spectral gap and the spin Chern numbers are well-defined and non-trivial. These conditions 
on being in the quantum spin Hall phase do not require time-reversal symmetry. The result 
materializes the general philosophy that topological insulators are topologically non-trivial bulk 
systems with edge currents that are protected against Anderson localization. 

1 Resume and discussion 

A quantum spin Hall system is a solid state phase of independent particles with spin. It has non- 
trivial bulk topology going along with edge states which carry spin currents and are not exposed to 
Anderson localization by a disordered perturbation. The first such system was suggested by Kane 
and Mele jKMTl IKM2] . These systems are the prime examples of so-called topological insulators 
[SRFL| IHK] . Until recently it was firmly believed (see all the reference just cited) that the quantum 
spin Hall phase is tied to a non-trivial Z2-invariant denned only for time reversal invariant systems 
with odd spin, but several newer works |SP} IYXS} ILSX] indicate that rather non-trivial spin Chern 
numbers |SWSHt IPro] are at the origin of the phenomena. In particular, time-reversal invariance is 
not a crucial requisite for the quantum spin Hall phase. Other than in the quantum Hall effect, there 
does not seem to be any physical bulk quantity directly connected either Z 2 -invariants or spin Chern 
numbers. However, it is usually argued that non-trivial invariants of the bulk lead to non-trivial spin 
edge currents |SRFI4 IKM2} IHK] . A first rigorous proof of such a connection for periodic systems 
was given in our prior work |ASVj . Independently, Graf and Porta provide an alternative proof 
applying to more general situations |GPj . While interesting, this is definitely not yet satisfactory 
because disorder is an essential ingredient to quantum spin Hall systems and topological insulators 
in general. The main result of this work examines the stability of the edge currents under disordered 
perturbations: 
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Theorem For a disordered two-dimensional tight-binding model with spectral gap and approximate 
spin conservation, the spin edge currents do not vanish provided the spin Chern numbers are non- 
trivial. 

The concrete model, explicit hypothesis and detailed result are described below. The argument 
is of perturbative nature around ^-invariant Hamiltonians (here s z is the z-component of the spin 
operator). Such Hamiltonians can be decomposed into a direct sum of Hamiltonians acting on the 
eigenspaces of s z . In presence of a gap, the Fermi projection of each such operator has a well-defined 
Chern number which is then called the spin Chern number of the total Hamiltonian. These Chern 
numbers lead to quantized edge currents by the argument for disordered quantum Hall systems 
[KRSj IEG] . If the spin- weighted sum of these currents does not vanish, it only remains to control by 
perturbation theory the effect of sMnvariance breaking terms. This is done in a quantitative manner 
in this paper. This may seem as a small addendum to prior works which, moreover, only faces minor 
technical difficulties, but it is nevertheless important because it is the first rigorous result supporting 
the conception of topological insulators widely accepted in the physics community. The argument 
also shows that vanishing of spin edge currents should rather be seen as a coincidence in the quantum 
spin Hall phase, and this even if disorder is larger. The result also agrees with rigorous findings on 
quasi-one-dimensional toy models for the edge states in a disordered environment [SSJ. 

In view of the above result and the argument leading to it, it is appealing to characterize the 
quantum spin Hall phase by non-trivial spin Chern numbers (which are not restricted to time- 
reversal invariant systems). As pointed out by Prodan [Pro], such spin Chern numbers can be 
defined whenever there is a gap in the spectrum of Ps z P (here P is the Fermi projection). Actually, 
this gap may remain open even though the gap of the Hamiltonian closes and becomes a mobility gap. 
In such a mobility gap, it is nevertheless possible to have well-defined edge currents if an adequate 
procedure for their definition is used (see |EGSj and the comment at the end of the paper). It can 
thus be expected that the above theorem also holds if the spectral gap condition is replaced by a 
mobility gap condition. In systems of physical interest there seems to be a true gap though and the 
role of the mobility gap is not as central as in the quantum Hall effect where no quantization would 
be possible without it. In view of all this, the spectral gap of Ps z P combined with non-trivial spin 
Chern numbers are the most adequate indicator of the quantum spin Hall phase. It is therefore of 
crucial importance to have efficient means to calculate the spin Chern number of a given model. For 
periodic sMnvariant systems the spin Chern numbers can be calculated by an efficient algorithm 
using only the finite size transfer matrices |ASVj . This is briefly reviewed at the end of Section [3J 

Acknowledgements: The author thanks Julio Cesar Avila and Carlos Villegas-Blas for the collab- 
oration |ASVj which led to the natural question examined in this paper, and Gian-Michele Graf for 
several interesting discussions on quantum Hall systems over the years, and topological insulators 
more recently. This work was supported by the DFG. 
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2 Models and hypothesis 



Let us consider a two-dimensional system of independent tight-binding particles with spin s 6 ^ 
on the square lattice Z 2 where over each point of Z 2 there are i?6l internal degrees of freedom. 
Hence the Hilbert space is £ 2 (Z 2 ) <g> © C r where r = 2s + 1, but we also set L = Rr and write 
C L = C R © C r . This setting allows to study particles on square, triangular and honeycomb lattice 
(see |ASVj . although this can probably also be found elsewhere). The Hamiltonian is a strongly 
continuous family (H w ) we Q of self-adjoint finite difference operators on £ 2 (Z 2 ) © C L indexed by a 
variable u, called the disorder configuration, taken from a compact dynamical system (O, T, Z 2 ) and 
satisfying the covariance relation 

H TaU = UaH^u:, aeZ 2 , 

where U a are the magnetic translations defined 

(U a (f>) n = e^ BnAa (j) n _ a , n A a = n x a 2 - n 2 a 1 . 

We also simply write H = (H^ujen (and this is then identified with an element of the C*-algebra 
of homogeneous observables used in |BESl IKRS] . but this will not be used here). The Hamiltonians 
H u are supposed to have matrix (hopping) elements of at most uniformly bounded finite range. The 
model is completed by fixing an invariant and ergodic probability P on Q. The main assumptions of 
the paper are the following: 

Gap hypothesis: There is an energy E g lying in a gap of the almost-sure spectrum of H = (H^uen- 

Approximate spin conservation: The norm of the commutator [H^,s z ] is uniformly bounded in 
norm by some sufficiently small constant C s < 1 smaller than the size of the gap. 

One consequence of these hypothesis is that there exists a homotopy H(X) = (H^X))^^ of 
covariant Hamiltonians between H(l) = H and H(0) commuting with s z such that the gap E g (X) 
associated to E g remains open along the homotopy (as C s is sufficiently small). It can be chosen as 

H{X) = H + l —±(s z Hs z -H) = H + ^[s z ,H]s z , X e [0, 1] . (1) 

By construction one has [H(0), s z ] = so that these two self-adjoint operators can be simultaneously 
diagonalized so that 

H(0) = F_ s (0) © ... © H„{0) , (2) 

where the direct sum decomposition is w.r.t. the eigenspaces of s z , that is £ 2 (Z 2 ) © C L = @i = _ s /hii 
where Hi = £ 2 (Z 2 ) © C^. This means that the Hamiltonians decomposes into r = 2s + 1 separate 
Hamiltonians, each of which can be studied as a quantum Hall system along the lines of [BESj IKRS] . 
This will be heavily used further below. 

Let us provide a concrete situation of physical interest where both the above hypothesis are 
satisfied, namely the Kane-Mele model on a honeycomb lattice with spin-orbit interactions. In this 
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situation, there are two internal degrees of freedom (R = 2) resulting from the bipartite structure of 
the lattice. The Hamiltonian on £ 2 (Z 2 ) <g> C 2r is now of the form 



H u = H hc + \ so H so + A Ra # Ra + A dis ^ y n \ n )( n 



neZ 2 



Here H bc is the discrete Laplacian on the honeycomb lattice, H so is the second nearest neighbor spin 
orbit coupled hopping, H Ks3 _ is the next nearest neighbor Rashba spin orbit hopping (all possibly with 
magnetic field), and finally the V n are independent random 2r x 2r matrices which are supposed to 
be real and commute with s z (typically they are diagonal with two random entries describing random 
potential values on two neighboring sites of the honeycomb lattice). A detailed description of the 
Hamiltonian is given in |KM2[ IAS V] . The main features of importance here are that [H so ,s z ] = 
and [iy Ra , s z ] 7^ so that the Rashba term is the only one breaking conservation of the s z component. 
It is known jKM2j that for s = ~, A Ra = A dis = and A so > the energy E g = lies in a gap of 
the spectrum and there is one band above and one below E g = 0. For sufficiently small A Ra and A dis 
(smaller, but possibly of the order of A so ) the above Gap hypothesis as well as the Approximate spin 
conservation both remain valid. Moreover, as discussed once again below, the model is a non-trivial 
topological insulator in this regime because the spin Chern numbers defined below do not vanish. 
Let us also point out that by multiplying A Ra by A G [0,1], one obtains another homotopy H U (X) 
from H u (l) = H w to H^(0) which then commutes with s z . This homotopy does not coincide with 
([rj, but may just as well be used below. 

3 Spin Chern numbers 

The spectrum of the spin operator s z is {— s, — s + 1, . . . , s — 1, s}, and each of these eigenvalues is in- 
finitely degenerate. Following Prodan [Pro], let us now consider the self-adjoint operator P(X)s z P(X) 
where P(X) = (P w (A)) w6 n denotes the spectral projection of H(\) = (H w (\)) ue n given in (JTJ) below 
the gap Eg(\). For A = 0, one has [P(0),s z ] = so that 



In particular, the spectrum of P(0)s z P(0) is still {— s, . . . , s} with each level infinitely degenerate. 
As A increases, this degeneracy is lifted leading to r = 2s + 1 spectral islands, but at least if C s 
is sufficiently small, there remain spectral gaps between these spectral islands. Hence it is possible 
to choose positively oriented contours 7; around each one of them and define a covariant family of 
orthogonal projections P/(A) = (^(A))^ by 



s z Pi(0) 



lP t (0) . 



(3) 





-i 



= —s 



s . 



Clearly, one has 



p(\) = p ^ ■ 



l=-s 
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As P U (X) has exponentially decreasing matrix elements due to the presence of the spectral gap, a 
standard argument (of Combes-Thomas type) shows that also P u> i(X) has exponentially decreasing 
matrix elements. Therefore the Chern number of P/(A) in the sense of |BESt IASS] is well-defined 
and an integer: 

Ch(P,(A)) = -^E P Tr £ (0|P^(A)[[^i,^(A)],[^2,P^(A)]]|0). 

As these integers are intrinsic to P(A) and its interplay with the spin operator s z , they are called the 
spin Chern numbers of P(A) and denoted by 

SCh / (P(A)) = Ch(P,(A)) , l = -s,...,s. 

By the additivity of Chern numbers, one has 

s 

Ch(P(A)) = SCh z (P(A)) . (4) 

l=-s 

Because A G [0, 1] h- >■ P(A) is a norm continuous path of projections, the spin Chern numbers are 
independent of A due to the homotopy invariance of the Chern numbers [BES] . For time-reversal in- 
variant systems, the Chern number Ch(P(A)) always vanishes because Ch(P(A)) = Ch(6P(A)0 _1 ) = 
— Ch(P(A)) if G is the time-reversal operator. This does not mean that all of the spin Chern numbers 
vanish though. It is worth mentioning that the spin Chern numbers are also well-defined if the gap 
assumption is replaced by a dynamical (Anderson) localization condition, but this is not used below 
where only edge states in the gap are studied. 

Of course, it is of crucial importance to be able to calculate the spin Chern numbers in a given 
model in order to determine whether or not one is in a topologically non-trivial quantum spin Hall 
phase. Due to the homotopy invariance properties of the spin Chern numbers discussed above, it 
is often sufficient to do this for a sMnvariant model without disorder. For simplicity let us also 
assume that the magnetic field vanishes (rational magnetic fields can be dealt with similarly, and 
small magnetic fields may again be turned of by a homotopy). Then the Hamiltonian decomposes 
into a direct sum ([2]) of periodic operators. If the hopping in the original model is nearest neighbor, 
then one such operator Hi acting on ^ 2 (Z 2 ) (g) is of the following form: 

H = nsi + t*s 1 + t 2 s* 2 + t;s 2 + t 3 s; + t;s 3 + w , s 3 = sts 2 . (5) 

Here Si and S 2 are the shift operators on £ 2 (Z 2 ) and the Tj, j = 1,2,3, and W = W* are Rx R 
matrices. If the hopping is not nearest neighbor, but finite range, this form can again be achieved by 
enlarging the number R of internal degrees of freedom (or in other words, enlarging the unit cell). As 
shown explicitly in [ASVj, operators on quadratic, triangular and hexagonal lattice can be brought 
into the form (jSJ). One now has to calculate the Chern number of the Fermi projection P; associated 
to Eg of one such operator Hi. The following procedure to achieve this is explained and justified in 
detail in |AS Vj . First build the family in k G [— tt, it] of 2R x 2P transfer matrices 

m f(E g l-(e lk T* + e~ lk T 1 + W)){T 2 + e lk T*)- 1 - (T 2 + e lk T*)* 
1 1 > \ (T 2 + e lk T*)~ l 
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under condition that the inverse of T 2 + e* fe T 3 * exists (this is required at least for almost all k). 
Because E g is in a gap of the spectrum of Hi, this transfer matrix is hyperbolic (no eigenvalues of 
modulus 1) and due to its symplectic symmetry it has exactly R eigenvalues inside of the unit disc and 
R reflected ones outside. Determine (numerically if necessary) R linearly independent generalized 
eigenvectors <p m (k) G <C 2R associated to the eigenvalues inside the unit disc and build an 2R x R 
matrix $(/c) = (</>i(fc), . . . , 0#(fc)). Then set 




where 1 and ±« are multiplied by the R x R identity matrix. The R x R matrix U(k) (and, in 
particular, the appearing operator inverse) is well-defined and turns out to be unitary. The associated 
winding number is equal to the desired Chern number: 

/ n dk 
— In det(U(k)) . 

This procedure can easily be implemented numerically [ASVj . For example, for the spin s — | 
Kane-Mele model with A Ra = A dis = one finds that the two spin Chern numbers of the lower band 
are 1 and —1. By homotopy, this remains true in the regime A so > max{A Ra , A dis }. 



4 Spin edge currents 

In order to study edge currents, we now simply restrict all operators to a half-space. Let II : £ 2 (Z 2 ) — > 
£ 2 (7j x N) denote the partial isometry from £ 2 (Z 2 ) onto £ 2 (Z x N), namely 1111* is the identity on 
£ 2 (ZxN) and 1TTI is the projection on f(ZxN) seen as subspace of f(Z 2 ). Then H(X) = {H w {X)) uen 
is defined by H W (X) = U.H U (X)H*. It is possible to include other covariant boundary conditions, but 
this leads to some further technicalities that are inessential. The family H(X) = (H^^q is still 
covariant in the 1-direction, but not the 2-direction. The observable for edge currents along the 
boundary is now the 1-component Ji(A) of the spin current operator which by Heisenberg's equation 
is given by Ji(A) = s z X\ = i[H(\), s z Xi\. Now let g(H) be a density matrix defined by spectral 
calculus from a real valued function g on R with J dE g(E) = 1. The edge current j c (g, A) is then 
defined by 

f(g,\) = f(Ji(A)p(F(A))) , (6) 

where T = T\ Tr 2 is the trace per unit volume in the 1-direction along the boundary and Tr 2 the 
usual trace in the 2-direction perpendicular to the boundary. For any operator family A = (A^^q 
on £ 2 (Z x N) which is homogeneous in the 1-direction, a formal definition of T is 

f{A) = E P Tr L (0,n 2 |X|0,n 2 ) . (7) 

n 2 >0 
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Further properties of T are collected in |KRSj (where T does not contain the trace Tr £ , but this 
leads to no changes). In particular, T is cyclic and the T-traceclass operators form a two-sided ideal 
in the algebra of covariant operators. As is obvious from this definition, not all covariant operators A 
are traceclass w.r.t. T due to the sum over n 2 in ([7]). Thus one has to prove that actually makes 
sense for sufficiently smooth functions g. This is part of the following main result of the paper. 

Theorem Let H = (H^)^^ satisfy the Gap hypothesis and the Approximate spin conservation. 
Then, if g G C 6 (M) is any positive function of unit integral supported in the spectral gap associated 
to Eg and P is the Fermi projection of H beneath this gap, 

s 

f(g,X) = ISCh(P) + XR(g,X), (8) 

l=~s 

where R(g,X) is a function satisfying for some constant C 

\R(g,X)\ < C\\g\\ 6 \\[s*,H}\\ . 

As already discussed above, the result shows that the spin edge currents are quantized for s z - 
invariant Hamiltonians and that they are non- vanishing if ^ s = _ s I SCh;(P) 7^ and both ||[s z , H}\\ 
and || (7 1| 6 are sufficiently small. In particular, the support of g cannot be chosen arbitrarily small for 
a given value of || [s z , H] \\ . 

5 Proof 

First of all, at A = the Hamiltonian decomposes into a direct sum (T5]). For each summand Hi(0) 
acting on £ 2 (Z 2 ) £g> C R , the Chern number of the Fermi projection Pi(0) below E g is Ch(Pi(0)) = 
SChi(P(0)). Therefore the results of [KRSl lEGj (which also contain the T-traceclass property of 
g(H(0)) that is shown once again below) imply 

f(t[H l (0),X 1 ]g{H l (0))) = SChKP(O)), 

where here T only contains a trace of C R and not C L . Due to ([3]), summing over / G {— s, . . . , s} 
shows that ([8]) holds for A = 0. Therefore, only remains to develop a controlled perturbation theory 
in A. 

Let us set H Q = H(0) and Hi = \[H, s z ]s z so that according to ([ID one has H = H + XH 1 with 
[H ,s z ] = and ||ifi|| < f||[-ff, s*]\\. The index u will be suppressed in the following. The first step 
will then be to decompose 

g(H(X)) = g(H(0)) + XG(X) , (9) 
where G(X) is some remainder. As </i(A) = Ji(0) + Xi[Hi, s z Xi], resembling the terms shows 

f(g^) = f(9,0) + xf{z[H 1 ,s z X 1 ]g(H )) + Af(J 1 (A)G(A)) . 
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The aim is thus to prove bounds on the last two summands. As in [EG], let us start from the 
Helffer-Sjostrand formula for an arbitrary self-adjoint and bounded operator H: 

g{H) = — I dxdy d^g(x,y) (z - H)' 1 , z = x + ty , (10) 

where, for some N > 1, 

n=0 

with some smooth, even, compactly supported function \ : (—1,1) — > [0,1] which is equal to 1 on 
[—8, 5]. The integral in ([TO]) is norm convergent and the function g, called a quasi-analytic extension 
of g, satisfies 

aa{x,y) = g( N+1 \x) X (y) + * f> (n) (*) ^ x'(y) , 

n=0 

so that, in particular, uniformly in x, y, 

\d,$(x,y)\ < C \\g\\ N+1 \y\ N . (11) 
We use ( TTU|) to calculate g(H(X)). Next let us replace the geometric resolvent identity 

' n — — t-— - n* + — L — (H(\)u* -uh(x)) — ^rrn*. (12) 



z-H(X) z-H(X) z-H(X) ' z-H{X) 

Replaced in (fTUjl . the first summand gives II g(H(X)) II* which vanishes because the support of g is 
contained in a gap of the spectrum of H(X). Thus the second summand leads to 

g(H(X)) = f dxdyckg(x,y) 1 (H (X) II* ~ U H (X)) JT7T\ ^* ■ (13) 

Let us first check that this positive operator is indeed T-traceclass by bounding its matrix elements: 
\(0,n 2 \g(H(X))\0,n 2 }\ < [ ^ \d*(x, y)\ |(0, n 2 \ L— (H(X) IF - U H(X)) ^r|0,n 2 )|. 

Now we need to control the decay of the matrix elements of the resolvents of H(X) and H(X). This 
can be done by a Combes-Thomas estimate which states that there is an r] > and C\ such that 

\(n\(z - H(X)y 1 \m}\ < e -v\v\\n-m\ _ 

Replacing this estimate and the corresponding one for the resolvent of H(X) shows that the matrix 
element |(0, ri2\g(H(X))\0, n 2 ) | is bounded above by a constant times 

I dxdy\dzg(x,y)\^ Yl e'^^ 2 -™ 2 ^ 71 ^™ 1 \{m\{H{X) IT - n H(X))\k)\ e -»2-fc 2 |-^]]*i| . 



mGZxN fceZ 2 



Now the operator H has only finite range hopping so that if (A) IT — UH(X) has non- vanishing 
matrix elements only in a finite distance away from the boundary (thus m 2 and k 2 are close to 0). 
Also m\ and k\ are only a finite distance apart. This shows 



\(0,n 2 \g(H(X))\0,n 2 )\ < C 2 fdxdy\d,g(x,y)\± £ £ 

n 2 >0 ^ R2 y n 2 >0 miGZ 



\n2\-r1\y\rn1 



Finally invoking the bound ( TTTT) with N = 4 and the fact that the integration domain is compact, 
one concludes that 

f(g(H(X)) < \(0,n 2 \g(H(X))\0,n 2 )\ < C 3 \\g\\ 5 . 

n 2 >0 



Next let us come to the perturbation theory in A. It is based on 

H(X) n* - n H(X) = (H U*-UH ) + \(HiTl* -XlHi) , 
and the resolvent identity 



(14) 



1 



1 



+ A 



1 



Hi 



1 



z - H(X) z-H Q ' " z-H z- H(X) ' 
as well as the resolvent identity for [z — i?(A)) -1 . Comparing with this leads to 



G(X) 



-1 

2^ 



dxdy d^g(x,y) 



+ 
+ 



z~H z- H(X) 
1 



(H(X) U*-UH(X)) 



1 



z-H 
1 



z-H, 



[Hi n* — n Hi) 

(H U*-UH 0/ 



z-H(X) 
1 



n* 



Hi 







z - H z- H(X) 



z-H(X) 



n* 



n* 



Because both g(H(X)) and g(H(0)) are T-traceclass due to the above, so is G(X). It follows that also 
Ji(X)G(X) is T-traceclass. Its trace will be bounded above again by a bound on its matrix elements. 
First of all, each of the summands contains Hi or Hi as a factor. Furthermore, each of the three 
summands contains one of the summands of (f!4"l) each of which is supported on a strip of finite width 
along the boundary. Combining this again with the Combes-Thomas estimates shows after similar 
estimates as above (with one supplementary factor 1 2/ 1 1 due to the extra resolvent) 

f(Ji(X)G(X)) < \(0,n 2 \Ji(X)G(X)\0,n 2 )\ < C 4 \\g\\ 6 \\Hi\\ . 

n 2 >0 

Hence it only remains to prove the bound 

T{i[Hi,s z Xi]g{%)) < C 5 \\g\\ e \\Hi\\ . 
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This follows again by proving an upper bound on the matrix elements (using the fact that Hi is only 
finite range so that the unboundedness of Xi plays no role). This finishes the proof of the theorem. 

Let us conclude with a few technical remarks. By pushing the expansions above a bit further, 
one can calculate the linear growth in (jSJ), namely the coefficient Ro(g) in R(g, A) = Ro(g) + 0(A). 
Actually, one can control errors to arbitrary order and calculate a formal expansion in A. Furthermore, 
in the situation where g(H(X)) ^ 0, the above calculation can still be carried through and edge 
currents can still be calculated by setting 

f (%(if(A))nV 1 (A)) = . (15) 

This is reasonable because for H(X) on £ 2 (Z 2 )£g>C L the current Ep((0\g(H(X)) Ji(A)|0)) in equilibrium 
indeed vanishes |BESj and thus (TTS1) just results by imposing that the expectation Ep is calculated 
before the sum over ri2 in the trace T. The remaining other term is then bounded by the above 
arguments. This procedure is similar to the one in [EGSJ. 
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